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ABSTRACT 


The work summarized in this report, which was carried out as a part of 
a NASA sponsored fissioning plasma research program, consisted of design 
power plant studies for two applications of the plasma core reactor: 

(1) As a breeder reactor 

(2) As a reactor able to transmute actinides effectively. 

In addition to the above applications the reactor produced electrical 
power with a high efficiency. 

A reactor subsystem was designed for each of the two applications. 
Tables 1 and 2 summarize the reactor design parameters for the breeder 
and the actinide transmuter, respectively. 

For the breeder reactor, neutronics calculations were carried out 
for a U-“233 plasma core with a molten salt breeding blanket. The 
primary objectives of the overall nuclear design were to design a 
reactor with a low critical mass (less than a few hundred kilograms 
U-233) and also a breeding ratio of 1.01. The later objective was a 
safety precaution to guard against diversion of fissionable material 
during blanket reprocessing. Since only enough U-233 would be bred in 
the blanket to replenish the amount depleted in the core, any diversion 
of U-233 during reprocessing would result in an insufficient amount of 
fissionable material to replenish the core and the reactor would shut 
down. Both of the above objectives were met in the final design. It 
is also possible to design for much higher breeding ratios in the range 
1 . 1 - 1 . 2 . 

The Plasma Core Actinide Transmutation reactor was designed to trans- 
mute the nuclear waste from conventional LNR^s. Each LWR is loaded with 
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Table 1. Plasma Core Breeder Reactor Reference Design 


Dimensions of Reactor Regions 


Plasma 

Helium 

BeO Moderator 
* 

Molten Salt 
BeO Reflector 


- 165 cm O.D. 

- 285 cm O.D. 

- 325 cm O.D. 

- 355 cm O.D. 

- 375 cm O.D. 


Fe Pressure Shell - 415 cm O.D. 


Critical Mass 
Breeding Ratio 
Power 


- 26.3 Kg 

- 1.0099 

- 2000 MWt 


Average Thermal Flux in Plasma - 3.42x10^^ neutrons 


cm^ sec 


Reactor Pressure 
Average Temperatures 
Plasma 
Helium 
Molten Salt 


“ 200 atm 


25,000 K 

- 3,000°K 

- 1,015°K 


Molten Salt Mass Flow Rate 


- 542 Kg/ sec 


Molten Salt Composition -71,7% LxF (99.995% Li^) , 16% BeF2» 12,3% TI1F4 


X 



Table 2. Reactor Characteristics of 

Plasma Core Actinide Transmutation Reactor 


Reference Design: 

Geometry: ‘ Spherical 

Dimensions of reactor regions: 


U233 plasma 

200 cm 

thickness 

He 

120 cm 

thickness 

Be moderator 

17 cm 

thickness 

Act. Oxide + Zr + He 

0.85 cm 

thickness 

Be reflector 

80-90 cm 

thickness 

Critical mass = 380 Kg 

Mass of actinides =1.27 metric tonne 


Power = 2000 MWt 

Avg. thermal flux in plasma 

= 2.06x1015 

n/cm^-sec 

Avg. thermal flux in actinides = 1.23x10 

14/2 

n/cm^-sec 

Reactor pressure = 200 atm. 

Temp: 


* 

^233 piQsjjja 

25000°K 


He 

3000°K 


Be moderator 

lOOO^K 


Act. Oxide + Zr + He 

800°K 


Be Reflector 

400-600°K. 



Actinide Composition: 74 atomic% Np^37. 7 atomic% Am^^l. 14 atomic% 

4 atomic% 
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88 metric tonne of uranium (3,3% and operated until a burnup of 

33,000 MWD/MTU is reached. The fuel is discharged from the reactor and 
cooled for 160 days. Next, the spent fuel is reprocessed during which 
100% of Np, Am, Cm, and higher actinides are separated from the other 
components. The concentrations of these actinides are calculated by 
ORIGEN and tabulated. These actinides are then manufactured as oxides 
into zirconium clad fuel rods and charged as fuel assemblies in the 
reflector region of the plasma core actinide transmutation reactor. 
Results of actinide burnup calculations for an equilibrium plasma core 
transmuter servicing 27 PWR^s show that after 13 cycles the actinide 
inventory has stabilized to about 2.6 times its initial loading. There 
are two mechanisms for the removal of actinides: 

(1) They are fissioned directly in the plasma core actinide 
transmuter 

(2) They are removed as U or Pu. 

The U and Pu can be used in other reactors. In the equilibrium cycle, 
about 7% of the actinides are directly fissioned away, while about 31% 
is removed by reprocessing. 

Fluid mechanics, heat transfer, and mechanical design considerations 
for both reactors are described in the report. 

Since it is desirable to have the Plasma Core Breeder Reactor (PCBR) 
be a self-contained unit, generating its own new fuel, an on-line repro- 
cessing system for the molten salt blanket is a necessity. Chapter 8 
describes protactinium removal and salt purification processes, calcula- 
tions of expected flow rates, and equilibrium concentrations of various 
isotopes present in the system. 



In order to achieve maximum effectiveness from the high temperature 
coolants from either of the two plasma core reactors, it was decided 
that a ternary power cycle would produce the highest efficiency power 
plant • The ternary cycle consists of a combination of MHD, gas turbine, 
and Rankine cycle energy conversion units. Two concepts were investi- 
gated systems with and without a high temperature regenerator in the 
helium loop. 

The achieved objectives of the study were as follows: 

(1) Model the nuclear IlHD power plant cycle. 

(2) Analyze the power output from the three energy conversion 
units and evaluate plant overall efficiency. 

(3) Make a parametric study of the effect of changing operating 
variables on plant overall performance. 

All studies used values for Input data according to current commercial 
technology (i.e. efficiencies for steam cycle components, gas turbine, 
and compressors) or with current use in MHD research. 

The modeling of the MHD cycle consisted of defining a pseudo- 
Bray ton cycle and treating the expansion within the MHD generator in 
a similar manner as in a gas turbine. In order to analyse the two 
systems it was necessary to write two computer codes: 

(1) HMHD-1 code to analyze the nuclear MHD power plant without 

regeneration in the helium loop 

(2) NMHD-2 code to analyze the nuclear MHD power plant with 

regeneration in the helium loop. 

Table 3 lists input parameters for each system. 

A study was made of the effect on overall efficiency of varying the 
reactor coolant outlet temperature from 3000^K to 4000\ for the two 
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systems* Table 4 lists typical results, showing an overall plant ef- 
ficiency as large as 70%. 

For Nuclear MHD Power Plant with regeneration, the major contribution 
of the electric power is produced in the top of the power cycle by the 
MHD subsystem (33.97% - 45.49% from 100% heat produced by the reactor). 

The power productxon has been shxfted toward the top of the ternary cycle 
with a positive effect on overall efficiency. This system produces 
overall efficiencies that are 25-35% higher than actual power plants in 
use and that are 15-20% higher than the expected coal-fired MHD power 
plants. 

For Nuclear MHD Power Plants without regeneration, the major contri- 
bution of electric power is due to the steam turbine subsystem (36.03% - 
36.36% from 100% heat produced by the reactor). Due to a significant 
fraction of the electric power being produced by the steam cycle with a 
low efficiency (40%) , it is desirable to shift the power production 
toward the top of the cycle to improve the overall efficiency. This can 
be achieved by reducing the mass flow rate of helium within the inner 
loop and increasing the pressure ratio of the MHD generator. This system 
produced overall efficiencies that are 15-20% higher than actual power 
plants in use and that are 5-10% higher than the expected coal-fired MHD 
power plant. Due to the relatively low temperatures within the helium 
loop, this type of power plant could be considered as a first step in a 
national program of implementation of MHD poxfer plants with a nuclear 
source. 
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Table 4» Plant Net Overall Efficiencies For MHD Inlet Temperature Variation. 
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1. INTRODUCTION 


As part of its policy of supporting research and development programs 

which reside on the frontier of power technology, the National Aeronautics 

and Space Administration has sponsored work in gaseous fueled reactors and 

plasma research. The original thrust of the NASA sponsored research, aimed 

toward development of a space propulsion engine, led to two gas-core reactor 

concepts - the light bulb and the coaxial flow nuclear reactor concepts.^ 

Although budgetary and policy factors terminated the development of nuclear 

powered propulsion engines, the concept of a UF^ fueled gas core reactor was 

shown to be very attractive for several other applications » 

NASA has continued supporting an ongoing fissioning plasma research 

program consisting of cavity reactor criticality tests, fluid mechanics 

tests, investigations of uranium optical emission spectra, radiant heat 

transfer, power plant studies, and related theoretical work.^^ These 

studies have shoxm that UFg fueled reactors can be quite versatile with 

respect to power, pressure, operating temperature, and the modes of power 

extraction. Possible power conversion systems include Brayton cycles, 

12 13 14 17 18 

Rankine cycles, MHD generators, and thermionic diodes. > ^ > Recent 

results of research on the pumping of lasers by fission fragment inter- 
actions w^ith a laser gas mixture indicate the possibility of the power 

10 12 

extraction in the form of coherent light, * Another potential applica- 
tion of the gas core reactor is its use for nuclear waste disposal by nuclear 

10 12 19 20 

transmutation (Gas Core Actinide Transmutation Reactor, GCATR) , ’ ^ * 

Recent work sponsored by NASA at Georgia Tech on the Gas- Core 
Breeder Reactor was reported in References 17 and 18. Further work on 
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the Gas Core Actinide Transmutation Reactor was reported in Refs. 19 
and 20. 

This semi-annual report summarizes results of work performed from 
March 1, 1977 to August 31, 1977 and NASA Research Grant NSG-1288. Work 
was performed in connection with the UF^ fuel under near "state-of-the 
art" temperature conditions, and also on the high temperature fissioning- 
plasma co-axial flow scheme. This report contains results for the application 
of the high temperature fissioning-plasma core to transmutation and 
breeding. 

Chapters 2 to 5 apply to the plasma core transmutation reactor and 
Chapters 6 and 7 relate to the breeder. Chapter 8 applies to both the 
MHD generator, a component in both systems, and Chapter 10 encompasses 
the system designs for both applications, showing why the fissioning 
plasma system is so extraordinarily attractive. 
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1 . NUCLEAR ANALYSIS OE THE PLASMA CORE ACTINIDE TRANSMUTATION REACTOR 


The objectives of the nuclear analysis of the plasma core actinide’ 
transmutation reactor are: 

(1) design a reactor system capable of producing power; 

(2) design a reactor able to effectively transmute actinides. 

The spent fuel discharged from a LWR consists of structural materials, 
unfissioned uranium, converted plutonium, other actinides, and fission 
products. The ratio of these components by weight is as follows: 

structural : uranium : plutonium : fission products : other actinides 

256 1023 9 36 1 

Despite the fact that the other actinides is the smallest component,, they 
are very long lived. After 10^ years, most of the other components will 
have decayed to stable isotopes, but these actinides will still be radio- 
active and may be a significant health hazard in the future. This is the 
rationale for putting these actinides in a reactor to transmute them to 
short lived fission products. 

The transmutation strategy used for the present calculations is 
shown in Fig. 2.1. The analysis was performed using the cross section 
code^ MG^, the multi-group diffusion code^ MACH-I, and the isotope deple- 
tion code^ ORIGEN. The flow diagram for the analysis is shown in Fig. 2.2. 

Each LWR is loaded with 88 metric tonne of uranium (3.3% U^^^) and 
operated until a burnup of 33,000 MWD/MTU is reached. The fuel is dis- 
charged from the reactor and cooled for 160 days. Next, the spent fuel 
is reprocessed during which 100% of Np, Am, Cm, and higher actinides are 
separated from the other components. The concentrations of these 
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Fig. 2.1 Strategy for actinide transmutation 
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actinides are calculated by ORIGEN^ and are tabulated in Table 2.1. 

These actinides are then manufactured into fuel rods and charged into 
the plasma core actinide transmutation reactor. 

The basic core configuration is shown in Fig. 2.3. 

For simplicity, spherical geometry is used. There are five regions: 

(1) Region I consists of the hot fissioning uranium 233 plasma. 

The bulk temperature is assumed to be around 25,000\ at a pressure of 
200 atmospheres. Because of the ionization of the uranium atoms at such 
temperatures, the plasma density is lower than that predicted by 

the perfect gas law. The equation of state of uranium at such tempera- 
tures and pressures is given by Ragsdale.^ 

(2) Region II consists of a helium layer# The bulk temperature is 
assumed to be at 3,000 K at 200 atm# The perfect gas law is assumed to 
be valid for helium at these conditions. 

(3) Region III consists of a solid liner at 1000 K# For a fast 
system, stainless steel is chosen as the liner material. For a thermal 
system, beryllium is used to act as reflector and ^moderator. 

(4) Region IV consists of He cooled, Zr clad actinide fuel rods at 
800 K, and 200 atmospheres. The actinides are assumed to be present as 
oxides. Only the principal actinides, Np-237, Ain-241^ Am^243 and Cm-244 
are included. The other actinides are very small. The concentrations 

of this region by volume is assumed to be 43% actinide rods, 12% Zr clad- 
ding, and 45% He coolant. 

(5) Region V consists of a reflector* For a fast system, iron is 
used; for a thermal system, beryllium is chosen. 
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Nuclide 

Gran^-atoms 

Kg. 

Np237 

-2;04+4 

4.82+3 

Np239 

3.27-3 

7.81-4 

Am241 

1.93+3 

4.65+2 

Am242m 

3.87+1 

9.36+0 

Am242 

4,65-4' 

1.12-4 

Am243 

3.09+3 

9.44+2 

Cm242 

4.50+1 

1.09+1 

Cm243 

3.22+0 , 

7.82-1 

Cm244 

1.17+3 

2.86+2 

Cm245 

7.89+1 

1.93+1 

Cm246 

9.00+0 

2.22+0 

On247 

1.16-1 

2.86-2 

Cm248 

7..7B-3 

1.93-3 - 

Cm250 

4.49-11 

1.12-11 

Cf249 

' 7.01-5 

1.74-5 

Cf250 

1.28-5 

3.20-6 

Cf251 

6.89-6 

1.73-6 

Total 

2.76+4 

6.57+3 


Table 2.1 

Actinide Concentrations Charged to Plasma Core Actinide Transmutation Reactor 


OEIGINAL 

OP POOR QUALITY 
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Both fast and thermal reactors were studied, but the reference design 
was chosen to be thermal in order to keep the critical mass at a reasonable 
value. In thermal reactors, beryllium is placed in region III and V. 

Table 2.2 summarizes the results obtained for the thermal plasma 
core reactor. 

Table 2.2. Critical Parameters Characterizing Thermal and PCATR 


Be thickness 

10 cm 

13 cm 

15 cm 

17 cm 

20 cm 

Critical mass 

22000kg 

2640kg 

562kg 

215kg 

53kg 

Critical radius 

778 cm 

383 cm 

229 cm 

166 . cm 

104 cm 

E (source) 

0.109 kev 

2.08 ev 

0.37 ev 

0.437 

0.35 ev 

«£a^> in act. reg. 

0.29 b 

0.245 h 

0.186 b 

0.222 

0.152 b 


He thickness - 25 cm.. Act. region thickness = 0.85 cm.. Outside Be 
reflector = 100 cm. 

It is observed that increasing the moderation available to core neutrons 
reduces the critical mass of the plasma core by a factor of 400. How- 
ever, as the core becomes more thermal, the average fission cross section 
in the actinide region also decreases. 

The dimensions of the reference plasma core actinide transmutation 
reactor is shown in Table 2.3. 

Table 2.3. Reference Plasma Core Actinide Transmutation Reactor 







Results of actinide burnup calculations for an equilibrium plasma 
core transmuter are shown in Table 2.4. It is observed that after 13 
cycles, the actinide inventory has stabilized to about 2.6 times its 
initial loading. There are two mechanisms for the removal of actinides 
(1) they are fissioned directly in the plasma core actinide transmuter 
and (2) they are removed as U or Pu. The XJ and Pu can be used in other 
reactors. In the equilibrium cycle, about 7% of the actinides are 
directly fissioned away, while about 31% is removed by reprocessing. 
This situation is illustrated in Table 2.5. 
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Table 2.4 


Actinide Burnup in Plasma Core Actinide Transmutation Reactor 

of irradiation, 265 days of cooling, 365 days of reprocessing 
^ Daughters) and fuel fabrication, 

s serviced (1.27 metric tonnes of actinides charged per cycle). 
THERM . 0.53648, EES = 1.035, FAST - 4.450, Avg. Thermal ElL -T.otih ' 


Batch 




Cycle No. | 

No. 

1 

2 

3 

4 

i_A_ 

' 6 

7 

8 

9 

10 

11 

12 

13 

1 

1.27 

.752 

.447 

.267 

1 

.160 

.099 

.065 

,045 

.033 

.026 

.022 

.020 

.018 

2 


1.27 

.752 

.447 

.267 

.160 

.099 

.065 

.045 

.033 

.026 

.022 

.020 

3 



1.27 

.752 

.447 

.267 

.160 

.099 

.065 

.045 

.033 

.026 

.022 

4 




1.27 

.752 

.447 

.267 ' 

.160 

.099 

.065 

.045 

.033 

.026 

5 





1.27 

.752 

.447 

.267 

.160 

.099 

.065 

.033 

.026 

6 






1.27 

.752 

.447 

.267 

.160 

.099 

,065 

.033 

7 


1 





1.27 

.752 

.447 

.267 

.160 

.099 

.065 

8 


i 






1,27 

.752 

.447 

.267 

.160 

.099 

9 


! 







1,27 

' .752 

.447 

.267 

.160 

10 










1.27 

.752 

.447 

: .267 

11 











1.27 

,752 

-447 

12 












1.27 

.752 

13 

* 

• 













1.27 

Total 

1.27 

2.02 

2.47 

2.74 

2.90 

3.00 

3.06 

TTii 1 

3.14 

3.16 

3.19 

3.20 

3.22 


Table 2*5 Actinide Xnveiitory During Equilibrium PCATR Cycle 


Beginning 
of Cycle 

3.3MT 

Fissioned 

- 0.23MT (7%) 

Reproc . 

“ 1.02MT (31%) 

End of 
Cycle 

2.05MT 

Charge 

1.27MT 

Beginning of 
Next Cycle 

3.32MT 
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3. FLUID MECHANICS OF PLASMA CAVITY REACTORS 


In plasma core reactors (PCR) the gaseous uranium is confined in 
a somewhat spherical geometry by the hydrodyn^ic forces exerted by 
the coolant. The coolant enters radially towards the central uranium 
plasma, providing both containment, and cooling. Due to the nature of 
processes inherent to PCRs, there are some very interesting fluid 
dynamic problems. The gaseous uranium metal is not completely con- 
fined by the coolant and, therefore, moves slowly through the cavity 
and is exhausted with the coolant. As the plasma moves through the 
core it radiates its energy to the coolant, which leaves the reactor in 
the neighborhood of 3000°K, 

Figure 3,1 shows a schematic of the flow expected to be encountered 
in a PCR, The advantages and the disadvantages of the PCR all stem from 
the gaseous state of uranium fuel. By being in the gaseous form the 
maximum operating temperature is increased by 10 fold over conventional 
power sources. Also, very small critical masses are possible. One of 
the main disadvantages is that the fuel moves through the reactor as 
shown in Fig, 3,1, The amount of fuel that goes through the system is 
very important. Small core loadings are of no value if the mass flow 
rate of uranium is excessively high, and if a large piping system is re- 
quired to be full of expensive fully enriched uranium. Therefore, one 
of the first goals of any PCR fluid dynamic analysis and design would be 
a small fuel to coolant flow ratio. 

In the late 1960^ s and early 1970’ s a considerable amount of work 
was done xd.th hydrodynamic containment schemes. These consisted of 
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rotating cylindrical flows, vortex flows, coaxial flows, and spherical- 
radial flows. Experiments and analysis were done on all flow arrange- 
ments to identify trends and establish parameters important to good 
confinement and heat transfer. The results of two of these efforts are 
shown in Fig. 3.2. The data shown are for both cylindrical and spherical 
coaxial flows. The mass flow ratio shown as a function 

of the ratio fuel volume/cavity volume. Both graphs show that for a 
spherical geometry a high mass flow ratio (above 100/1) is obtained only 
when the volume ratio is less than 0.25. This implies a radius ratio 

-i/t ^ T of less than 0.63. For the reactor core sizes studied in 
fuel' total 

this report the above finding leads to the following conclusions for 
cavity size: 


CORE I 

CORE 20 

Fuel 

Radius 

"'50 cm 

~200 cm 

Fuel + coolant 
radius for 

m /m = 
s e 100 

~95 cm 

"'320 cm 

ik - 

for ^ ~ ^ 

m 1000 

c 

"'109 cm 

~430 cm 


Reactors of the size listed above would be of about 2000 MW(th) and range 
from about y to 2 times the size of present day PWR pressure vessels. 

The exact nature of the flow pattern in the core is unknown at this 
time. Knowledge of this would require extensive experimentation and 
analytical work to solve the coupled energy and momentum equations. For 
this project it is felt that identification of the main design goal 
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Fig. 3.2(a) Porous Wall Data on Gas Core 
Performance Map (A) 
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Fig. 3.2(b) Comparison of Coaxial-Flow and 

Spherical Cavity Containment Results ( 5 ) 
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(low fuel flow rate) and the parameters which control the attainment 
of this goal is sufficient. The problems of core temperature profiles 
and wall cooling needs will be discussed in later sections. 

For the plasma core reactor the coolant will have to absorb the 
thermal radiation emitted by the uranium plasma and transfer this energy 
to a MHD generator, gas turbine, and various heat exchangers. Due to 
the high operating temperature of the uranium plasma ('"25,000^K) and the 
MHD cycle, the coolant used in the PCR will have to have some additional 
characteristics above and beyond those of conventional coolants. Table 3.1 
shows some of the properties which the PCR working fluid must have. Items 
6-7-8 are the "extra*' characteristics required by the PCR. 

When the PCR was being considered by NASA as a propulsion device, 

16 7 8 

hydrogen was used as the coolant because of its high specific impulse. ^ ^ ^ 

However, for central station power (on earth, other planets, or a space 

station) the dangers and chemical reactivity of hydrogen seem to make other 

choices more attractive. Also, hydrogen is optically transparent to radiation 
below 1216£ until a temperature of 5000 to 6000 is reached. As Fig. 3.3. 

shows, the emission spectra of uranium plasmas is in the range of 3000 to 
8000 jS , and hence, hydrogen requires an added submicron sized seed to 

increase its absorption. Unfortunately, most gases suffer from low absorp- 
tion in the range of interest and will also require a seed. 

Other coolants for central station power plants are listed in Table 3.2. 
Hydrogen and nitrogen can be discarded immediately because of their high 
chemical reactivity. Carbon dioxide is good from a thermodynamic and heat 
transfer point of view, but chemical reactivity and decomposition at high 
temperatures could cause numerous problems. Helium and other rare gases 
are good because of their chemical, nuclear, and temperature stability. 
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Table 3.1 Coolant Propertias for the PCR 


A. Physical and Nuclear Characteristics 

1) Low neutron absorption cross section 

2) Low Induced radioactivity 

3) Good radiation stability 

4) Good thermal stability 

5) Compatibility with structural and component materials at high 
temperature 

6) Large thermal radiation absorption cross section 

7) High electrical conductivity @ 3000°K 

8) Molecular or atomic weight much lower than uranium 
( to enhance separation) 

B. Economic Characteristics 

1) Resumable cost 

2) Good availability 

3) Low pumping or compressing power requirements 

4) High thermal conductivity and small viscosities 


Table 3.2 Possible Coolants for the PCR 

1) Hydrogen 

2) Helium 

3) Carbon Dioxide 

4) Nitrogen 
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REtATIVE INTENSmr 



Fig, 3,3 Spectral Distribution of Radiation 

from Various Uranium Light Sources (10) 
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These gases also have acceptable thermodynamic and heat transfer pro- 
perties. Helium, due to its extensive operational experience, low 
atomic weight and availability is the final choice. 

Helium has been shown to readily meet, the requirements set forth in 

1 9 11 

Table 3.1 , except for items 6 and 7. ’ ^ Helium, like other gases,, 
has line absorption being the predominate absorption mechanism; thereby 
only absorbing a substantial amount of radiation at certain wavelengths. 
Since the width of these lines is small, the integrated absorption over 
a, range of photon energies would be negligible. Exact data for helium 
absorption in the range emitted by the PCR has not been found. A litera- 
ture search revealed data only do^ to lOeV ('"124£) . Figure 4 shows data 
in which helium absorption is compared to that of hydrogen (for which 
data in the PCR range of interests is available) . The data shows that the 
helium absorption is a little higher than hydrogen in the low energy 
ranges. Evidence of seeded helium radiation absorption experiments were 

found, but energy content of the seeded aerosols was reported instead of 

13 

absorption cross sections. Figure 3.5 shows typical results of these 
experiments. Here the helium seeded aerosols tend to show a somewhat 
smaller absorption than hydrogen. 

Since no data were found for helium absorption, a few assumptions 
were made which enabled the study to continue. 

1) Since the Lyman Series for helium ranges from 230 to 30o£, it 
was assumed that discrete absorption would take place in radiation 
fields of longer wavelength, therefore, necessitating seeding no matter 
what the absorption coefficient. 

2) Since all data found indicated helium absorption 
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Fig. 3.4 Photon Cross Sections for He 
and Hydrogen (12) 
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of low energy radiation to be *' close” to that of hydrogen, and since 
the absorption in a seeded gas is governed by the seed more so than the 
gas, data for absorption coefficients in hot seeded hydrogen will be 
used for this report. 

The purpose of seeding the helium is twofold: to increase its 

absorption coefficient and to make this coefficient as independent of 
wavelength as possible. Therefore, the first requirement for a seed 
material is a high, wavelength independent absorption cross section. 
Secondly, the seed material cannot readily react with the coolant, uranium, 
or any structural material. Also, the seed should not agglomerate as 
this causes a decrease in uniformity and decreases the absorption ef- 
ficiency. These requirements are listed in Table 3.3. 

A large amount of work, both theoretical and experimental, has been 
carried out with seeds of carbon, tungsten, iron, and silicon. Carbon 
was originally dismissed because of its reactivity with hydrogen, but 
with helium as the carrier gas this should not be a problem. Figure 3.6 
shows a comparison of theoretical absorption coefficients for tungsten, 
silicon, and carbon at 2000X as a function of particle size. Figure 3.7 
shows the attenuation coefficient of a hydrogen— carbon aerosol at 
3450°F. 

Table 3.3: Seed Material Requirements 

1. Good absorption; independent of wavelength 

2. Chemically non-reactive with PCR materials 

3. Does not agglomerate easily 

4. Compatible with MHD power generation 

5. Easily introduced into the helium gas. 
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ABSORPTION COEFFICIENT <>'%n 



Fig. 3.6 Theoretical Absorption Coefficient of Tungsten, Silicon 
and Carbon as a Function of Particle Radius 
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MASS ATTENUATION COEFFICIENT 



WAVELENGTH- ANGSTROMS x 10*® 


. Fig. 3.7 Attenuation Coefficient of Carbon-Hydrogen Aerosol at 3450®F 
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Comparing this with Fig- 3,8 , which shows the attenuation coef- 
ficient for tungsten aerosols at 3000^F, one can see that carbon aerosols 
have a clear advantage. Also, Fig. 3.7 shows that the carbon seeds 
produce the required independence with wavelength and have a high over- 
all absorption coefficient, being around 5*10^ cm^/gm. For these reasons 
carbon was chosen as the seed material for use in the PCR. 

The overall ability of the gas-seed mixture to attenuate thermal 
radiation is a function of not only the absorption coefficient of the seed, 
but also the density of the seed. The density of the seed is in turn 
limited by the aerosol generator capabilities, the degree of agglomeration, 
and the particle size. 

Aerosol particle densities from 4x10“^ to 8x10*"^ gm/cc have been suc- 
cessfully produced. From this range of particle densities the attenuation 

1 ic 

parameter RCcm^-*-) can be calculated for various seed-to-coolant density 
ratios. The results of these calculations are presented in Table 3.4. 

The amount of seed needed in the PCR will be determined by the heat 
loading on the liner wall. Since the core liner will be composed of 
moderator and structural material, the temperature and heat flux limi- 
tations on this component will have to be compatible with the material 
properties. Beryllium is proposed as the moderating material and stainless 
steel is for the structural material. Use of these materials will limit 
the temperature and heat flux to values somewhat close to those of 
present LMFBR designs (600^K and 2.5 . Figure 3.9 shows the fraction 

of one core power deposited in the liner. For the two core configurations 
_ 

R is the parameter for exponential attenuation of radiation in the formula 
l/r^ - where I is the Intensity, 
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Fig. 3.9 Fraction of Power Deposited in the Core Liner 
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studied, a value of R from 0.1 to 0.5cm~^ would result in a very small 

heat load on the liner. Table 3.5 summarizes the thermal radiation 

absorption results for the two cores studied. The fraction of the heat 

deposited in the liner ranges from negligible to 13.5 MW with a maximum 

heat flux of 1.07 MW/m^. This value, although high, is completely within 

15 

today^s technology (the average heat flux in the CRBRP is '"1.27 
To obtain the attenuation shown in Table 3.5 would require particle den- 
sities in the range of 2*10~^gm/cc to !• lO^^gm/cc. These are well within 
current technological capabilities and the resulting particle to carrier 
densities are low-enough to keep agglomeration to a minimum. Operation 
in the geometrical and property ranges of Fig- 3.9 should solve the 
problem of wall heating in the core. The core exhaust nozzle and the 
MHI) duct and its associated piping will still require extensive calculations 
to determine the cooling required. 

Upon absorbing the radiant energy from the uranium plasma, the 
particle seeds will transfer this energy to the helium coolant. Unlike 
convent ial power sources which are limited in the possible AT which can 
be tolerated, the PGR^ since it has no structure in the core region proper, 
can produce any AT desired. The desired AT can be established by adjusting 
the mass flow rate of helium. The core inlet temperature is really the 
only fixed temperature in the entire reactor-power plant system. Since 
the core liner must be made of a material such as stainless steel and 
the actinide fuel rods must be clad in zircoloy, the core inlet temperature 
must be limited to somewhere around 600^K. The outlet temperature may 
now be adjusted to any desired level for efficient and reliable MHD 
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operation. Figure 3,10 shows the required mass flow rates for a given 
temperature rise. For efficient MHD operation the reactor exit temperature 
needs to be in the range of 2500°K to 3500°K, which implies a AT of 
~1900°K to 2900°K. For this AT a mass flow rate of 190 to 130 kg/sec 
would be required. For the sizes of the reactor cavity studied here 
this would result in incore velocities of 190 to 1600 cm/sec for a 50 cm 
radius core and 10 to 100 cm/sec for a 200 cm radius core. For efficient 
MHD operation an entrance Mach number of 0,5 is also required. For the 
MHD unit used in this design study, the entrance area is 0,0387m^, This 
implies an area ratio of (A/A*; where A* is the critical area size) of 1,524, 
Working backwards from to A^ (see Fig. 3.11) we can find that for a 

core exit Mach number of 0.1 (low Mach number is needed here to reduce 
friction and aerodynamic heating of the walls) we need an exit area of 
0.175m^ (or pipe diameter of 0.472nl). This will in turn give us a 
reactor exit velocity of 4. 20* 10** cm/sec. to 2. 65- 10** cm/sec. 

No attempt has been made to tackle the problem of the details of 
the flow stability and confinement of the plasma- coolant system. These 
would be out of the scope of this design project. Attempts have been made 
to ascertain an "order of magnitude" of the effects and parameters sped— 
cally associated with the FCR. Geometry, mass flow and coolant . property 
ranges which would be characteristic of a 2000 MW(th) central power 
statioUj. have been examined. A summary of the fluid mechanical aspects 
of the FCR design is given in Table 3.6. 
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erature Difference 



2800 3200 


Table 3 , 6 Fluid Mechanical Summar 


(Item) 

Coolant 

(Value) 
or name 

Helium 

Exit Temp Range 

3000 to 3500°K 

Fuel to coolant mass 
flow core 

1/100 to 1/1000 

Coolant thickness 

45 to '220 cm 

Plasma radius 

50 to 200 cm 

Coolant mass flow 

120 to 190 kg/sec 

Seed material 

carbon 

Seed size 

0 . 012 microns 

Seed density 

10“ ® to 10"^ gm/cc 

Coolant exit pipe 
diameter 

2. 65* 10“^ to 4.2-10*" 



















Reactor 

Exit 


MHD 

Inlet 



§ known: 

nil = 0.1 

= 3000 K 
calculated: 

= 0.175 


@ A 2 known: 

A 2 = 0.0387 
m 2 = 0.5 
calculated: 

A* = 0.0254 


Equation; 

A/A* = l/m[2/ it+iy,! +(jT-l)mV2J0 


Pig. 3.11 Area Calculation. Technique 
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4, HEAT TRANSFER IN THE PLASMA CORE REACTOR 


Plasma core reactors are capable of producing heat at extremely 

high temperatures for use in rocket propulson, MHD power generation, or 

process heat applications- Most of the work dealing with heat transfer 

in a gaseous core reactor has been concerned \n±th reactors used in 

rockets- Here the uranium plasma heated a hydrogen coolant that was 

used for propulsion* For central station power production the difference 

is that helium is used as a coolant and is exhausted through a MHD nozzie- 

Therefore, previous work on heat transfer of gaseous core reactors can 

be readily applied to this design study * 

The simplest case to analyze is when there is no mixing of the fuel 

1 

and coolant in the uranium-helium core region- In reality there will be 

some mixing and some convection effects at the outside surface of the 

uranium plasma- Therefore, the case analyzed would be a '’first estimate" 

of the temperature distribution in the core- 

In this study the steady state temperature profile as a function of 

radius in a spherical geometry is analyzed. Figure 4.1 shows the geometry 

and pertinent data- The helium and uranium gases are assumed to be grey 

gases, which means the radiation absorption coefficient is independent of 

wavelength- The containment wall is also assumed to be grey so that wall 

emissivity and reflectivity are independent of wavelength. The approach 

2 

used is that proposed by Ragsdale and Kascak. In this method the volume 
heat generation term (q"’ ) is assumed radially dependent and the absorption 
parameters temperature dependent. 
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Fig. 4.1 Geometry and Data 
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If one assumes that the heat flow is basically a diffusion process 


the heat flux in the core may be expressed by 


2 


q" (r) 


-4 dCqT^*) 
3k(T) dr 


(4.1) 


where q" = heat flux 


cr - Stefan-Boltzmann constant 

T = temperature 

r = radius 

k = absorption coefficients 


At steady state the heat flux can be related to the volumetric 
' heat generation q"’ (r) by 


q"(r) 


-r 


(r) 4Trr^dr/47Tr^ 


(4.2) 


Figure 4.2 shows calculated fission densities for these different 
PCR radii. If the fission density, and hence the heat generation rate 
is fitted by a polynominal such as: 


i'" = X 


1=0 




(4.3) 


then integration of Eq. 4.2 yields 


q" (r) = X ^ 


i+1 


x=0 1+3 


(4.4) 
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Substituting into Eq. 4.1 and rearranging yields 


M = - 3k(T) ^ ^ 

“ 16oT3 


(4.5) 


Equation 4*5 can be solved by numerical techniques with k(T) being 
supplied via experiment results* 

To solve Eq* 4.5 one must also have the boundary conditions for 
the edge temperature of the reactor (1^) • Unfortunately, this temperature 
is not known, but because of wall stress and creep limitations the wall 
temperature (T^) is known. This temperature will in turn affect the 
core edge temperature. To get a relationship between and T^ let us 
first consider a brightness temperature T.^^* T^ is defined so that 
gives the radiated heat flux at the core edge* is defined by the 
expression: 


where 


q" = aF(T^ - t'*) (4.6) 

b w 



= emissivity of the wall. 


Thus the brightness temperature is found to be: 



(4.7) 
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i+1 


n q . r 

where q” = - edge heat flux 


x=o 


To determine the edge temperature from T, , Ragsdale and Kascak pro- 

D 

posed the following relationship: 


T 

e 



2r^k(T)j 


(4.8) 


Since k(T) is a function of T and, therefore, must be evaluated at 
Eq. 4.8 must be solved iteratively. The use of the secant method 
provides quick convergence. 

Once is known T(r) can be solved by any appropriate numerical 
method. A fourth-order Runga-Kutta program was written and implemented 

for this purpose. The data for k(T) was input into the program as a 

polynomial fit to the data shown in Rig. 4.3. The program was checked 

1 2 

for accuracy against previous work ’ and showed good agreement. 

Temperature distributions for the design basis core of 200 to 
20 cm in radius were calculated and are shown in Fig. 4.4. The com- 
bination of lower temperature (and hence higher absorption coefficient) 
and a higher fissioning rate at the core edge cause the very rapid 
temperature increase in the first few centimeters inside the core. This 
rapid rise could also mean that the core edge temperature has been under- 
estimated; however, the calculations were checked and the values used 

for T- and T were consistent. These calculations showed that maximum 
b e 
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Temperature (T/T 



Fig. 4.4 Core Temperature vs Radius 
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core temperatures ranged from 21,400 to 21,800'^K with an average 
temperature in the neighborhood of 18,500°K. 

The transuranium actinides (neptunium, plutonium, americium, 
curium, etc.) will be introduced into the PGR in the form of zirconiuri 
.alloy clad rods located just outside the beryllium reflector. The 
actinides will be in the form of oxides and will be cooled by the helium 
coolant used in the core. The design constraints for the actinide rods 
will have to be the same as those used for present day power reactors for 
safety reasons. These are: 

Linear heat rate. . . 500 watts/cm. . . to control amount of stored 

energy in the fuel. 

Maximum clad temperature. . . 500°C. . . to maintain clad strength. 

There is a scarcity of thermal-physical property data on the oxides 
of transplutonium actinides and there is no data on the behavior of 
mixtures. Their known melting points, are high (~2400*^C) and actinide 

4 

oxide densities are all about the same being "'ll gm/cc. Consequently, 
since little is known about the actinide oxides , and a sizable fraction 
of the actinide rods will be uranium oxide, it is assumed that their 

5 

physical properties are the same as UO^. 

A high burnup rate in the actinide' fuel rods is desirable so as 
to shorten the time required for transforming the actinides by fission. 
The burnup rate is directly proportional to the volumetric heat gene- 
ration rate, which is: 
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R. 

X 


^fisslol 


(4.9) 


(r) 


-t. 


where 

= heat generation rate, watts 

R^ = energy released per fission (Joules) of ith isotope 

^fission ” fission density of ith isotope, fissions /cm^-sec. 


The value of R^ for the actinides is not known, but it should be close 
to that of uranium or plutonium ('"200MeV/fission) • Consequently, the 
volumetric heat generation for the actinide rods is : 

MeV ^ ^ f 1 Joule W - fissions\ 

fxssxony \,6. 242x10^2 MevJv cm^ sec>’ 

For the initial actinide fuel loading into the PCR 

^fission “ 4 •34x10^^ fissions/cm^-sec 



which gives 

q'” = 13* 90 watts /cm^ . 


This value is considerably lower than the heat rate in present LWRs so 
the probability of exceeding one of the design limits is very small. 

The geometry for the actinide fuel rods was chosen to give a 
minimum actinide region thickness, as this helped give low critical masses. 
This also helped keep the flux flat in the actinide region (for even 
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burn-up), and helped to limit the volume of actinrdes required to fill 
the region (at an inner radius of 3,4 meters even a small thickness 
has quite a large volume) « The radius for the fuel rods was set at 
0.268 cm with a 0.035 cm thick clad and a gap thickness of 0.015 cm. 

This gave an overall radius of 0.318 cm (0.25 in. diameter). The rods 
will be spaced with a 0.142 cm diameter wire wrap with a pitch of 25 cm. 

The actinide fuel, rods are one meter long and positioned just out- 
side the beryllium liner as shown in Fig. 4.5. Calculations showed, that 
for the dimensions shown in Fig. 4.5, the neutron flux variation in 
going from the center of the actinide rods to either end is only 1.0%. 
Hence, it is a good assumption to assume that the flux and, thus, the 
power does not change over the rod length. Therefore, 

q"’ 7 ^f(r,z); q”^ = const - 13.90 watts/cm^. 


With this value we can now calculate the helium temperature at the 
exit to the actinide fuel rods. Unfortunately, this is an iterative pro- 
cedure because the design point temperatures are the core outlet and the 
actinide inlet. Thus the mass flow rate will have to be varied until 
the design core outlet temperature is met. The procedure will be: 


1 ) 

2 ) 

3) 

4) 


Assume mass flow rate mx 

Trr^ Hq”* ^n 

Calculate T^ = T “T, + f^ 

out actinxdes o In 


Calculate T == T 

out core o 


== T + 

""a 

a p 


Check if T - 3500 K if not, 

o 

c 

redo calculation with new mass flow m 2 . 
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The above procedure must be done within the design limits set forth 
and with the following data: 

P = 2000 m 
r^ = 0.269 cm 

H = 100 cm 

n = number of actinide fuel rods: 

= 2725 for first fuel loading 

T. = 642°K 

xn 

a 

^ — C OA 

^helium * gm^K 

q"^ = 13,9 watts/cm^ 

r = 0.318 cm 

w 

The calculation must also be done taking into account the effects of 
the uncertainties due to manufacturing tolerances, physical property 
variations and flow maldistributions. The effects these uncertainities 
have on the temperature can best be accounted for by applying hot channel/ 
hot spot factors to the computation.' Because of the similarity of the 
actinide rods and the fuel rods of the Clinch River Breeder Reactor the 

g 

same hot channel/hot spot factors used in its design will be used here. 

The nuclear hot spot factor for axial variation in the flux will be 
ignored here due to the arguments given earlier for assuming a flat 
heat' generation rate in the actinide fuel rods. The factors used are: 


Coolant 

Film 

Heat FIuk 

F 

Ah 

F 

AT 

F 

q 

1.232 

1.168 

1.081 
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Using these factors the actinide fuel rod heat transfer calculations 
yielded the results shown in Table 4.1. As can be seen the maximum 
clad temperature turns out to be 647^0 (374^0) and is still well within 
the design limit of SOO^C. These calculations are for startup of 
actinide transmutation. Tl^e calculated tem^-ratures at later cycles 
will probably be higher. 


Table 4.1 Actinide Region Thermal Parameters 


Maximum CLAD Temperature 
Inlet Helium Temperature 
Outlet Helium Temperature 
Mass Flow Rate 
Helium Velocity 


646°K 
642°K 
643. 7 

134.65 kg/sec 

5.9x10^^ cm/sec (M = 0.4) 


These calculations represent scoping type work and are intended to 
give "order of magnitude" values for temperatures in each region. A 
more exact temperature profile for the core region would require detailed 
information about the flow field Inside the uranium plasma and in the 
helium coolant. Items that were not considered, such as fuel-coolant 
mixing and the energy of the uranium swept out of the core, need to be 
analyzed. Also, detailed pressure loss calculations need to be performed 
for the actinide and core regions. This work, however, is enough to 
permit a fairly accurate description of the operating conditions inside 
a PCR. Table 4.2 summarizes design points relevant to the heat transfer 
analysis of the PCR. 


OF pool 
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Table 4.2 Heat Transfer Summary 


ITEM 


VALUE 


Actinide Region: 

Helium Inlet Temperature 642°K 

Maximum Actinide Temperature 652°K 
Maximum Clad Temperature 646°K 

Helium Exit Temperature 643. 7°K 


Core Region; 

Mass Flow Rate 
Helium Exit Temperature 
Maximum Fuel Temperature 
Exit Helium Velocity 


134.65 kg/sec 
3500°K 
~21,400°K 
2.85x10^ cm/sec 


ORIGINAL PAGE IS 
05! POOR QUALITY 
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5 . DESIGN CONSIDERATIONS QE PLASMA CAVITY REACTORS 

Figure 5,1 is a schematic of the design base reactor* The relevant 
dimensions, obtained from the fluid mechanical, nuclear analyses, and 
thermal-hydraulic analysis ^ are also shoxm. The reactor is almost 3.40 , 
meters in radius and 7.20 meters in length. Table 5.1 gives the region 
volumes and components weight. 

The beryllium liner, besides serving as the reactor moderator, is 
also the porous wall and flow director for the helium coolant entrance 
into the reactor cavity. Beryllium was chosen for the liner because of 
its exceptional moderating capabilities. The helium temperature in the 
beryllium liner ranges from 370^0 to 400^^0 which is approximately one- 
fourth its melting point, so the metal should be able to withstand such 
a thermal load. There will also be very little pressure difference be- 
tween the liner inside and outside so the overall stress on the liner 
should be small. The beryllium metal is encased in a jacket of zirconium 
alloy and bolted under a small compressive load. Since the thermal ex- 
pansion of beryllium is greater than zircaloy, at operating temperatures 
the beryllium should be in even more of a compressive load that should 
equalize any pressure load on the beryllium. Figure 5.2 shows a detail 
drawing of the liner. 

The actinide fuel rods presented the first real problem in the over- 
all system layout. Since spherical geometry is best from a neutronic 

W 

and hydrodynamic standpoint for the core, the cylindrical actinide fuel 
rod did not really ”fit" anyx^here. It was decided to place the rods just 
outside the beryllium reflector in a cylindrical annulus region. Each 
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Table 5*1 Volume and 



(I) Uranium plasma 

(II) Helium coolant 

(III) Beryllium liner: 
Zirconium Alloy: 

Helium coolant: 

IV) Actinide rods : 

Zirconium Alloy: 

Helium coolant and bond: 

(V) Control drums : 

Cd: 

Beryllium Shield: 
Zirconium Alloy: 

(VI) Helium coolant outside 

of Beryllium liner 
and Actinide Region: 


*These are loadings for each yeat 


weights would be 3 to 4 times 


of Reactor Components 
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Zirconium alloy 
plate 
( . 635cm) 
typica 
coolant 

flow channels ^ 

(A .5 cm dia, hoi 
every 1.5 cm along 
outside) 



10 beryllium plates 


Typical tensioning 

bolts (1.27cm dia. bolt) 


Fig. 5.2 Liner Details 



■Additional actinide 
region plugged off for 
1st core loading, i^lugs 
^removed and additional 
rods inserted; at 
. equilibrium about 3 
rings of rods will 
be in core 


.035-*{lfr- 


\ I 

l•f--142 
-^l|<^015 « I 


Beryiium liner - 


Fig. 5.3 Actinide Fuel Rods 


All dimensions in cm 


Drive gears 


Zirconi 

can 



Moderator 
filled or 
empty 

Toward 

reactor center 




100 cm. 


f - fe; 


i. 




1 




Cd strip 


/ 

/ 


2.58crax.05cm 

deep 


A 


Fig. 5.4 Control Drum Details 
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rod with its associated clad and wire wrap will require 0*9194 cm width* 
The rods will be loaded in overlapping annulus regions as more actinides 
are introduced into the core* The initial charge of actinides should be 
about 445 Kg of actinide-oxides. This represents the yearly discharge 
from about 20 PWRs and requires one ring of actinide rods around the 
core. The equilibrium actinide loading would be about three rings of 
actinide rods or about 1500 Kg of actinide-oxides. Figure 5.3 shows the 
detail of the actinide region. 

The control drums will be zirconium alloy cylindrical cans filled 
with beryllium or beryllium oxide with a strip of cadmium on one side. 
Figure 5.4 shows a typical control drum. These drums will provide the 
reactivity control needed to assure adequate dynamic control during op- 
eration. 

A possible reactor layout is 'shown in Fig. 5.5. 
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6> NUCLEAR ANALYSIS-PLASMA CORE BREEDER REACTOR 


Neutronics calculations were carried out for a U-233 plasma core 
with a molten salt breeding blanket. The primary objectives of the over- 
all nuclear design were to design a reactor with a low critical mass 
(less than a few hundred kilograms U-233) and also a breeding ratio of 

1.01. The later objectives was a safety precaution to guard against di- 
version of fissionable material during blanket reprocessing. Since only 
enough U-233 would be bred in the blanket to replenish the amount deplet- 
ed in the core, any diversion of U-233 during reprocessing would result 
in an insufficient amount of fissionable material to replenish the core 
and the reactor would shut down. Both of the above objectives were met 
in the final design. 

The Mach-I code^ was used as the primary computational tool in the 
nuclear analysis. Mach-I is a one— dimensional diffusion theory code 

2 

which uses the 26-group ABBN cross section set of Bondarenko, et al. 

All neutronic calculations were performed by varying the plasma core 

radius to obtain K 1.000 with other dimensions held constant, 

eff 

Initially, four moderators were analyzed to determine the one most 
suitable for the final reactor design. The moderators selected were gra- 
phite, beryllium, beryllium-oxide, and heavy water. Properties of each 
of these are listed in Table 6.1. 

The geometry used for comparing the moderators is shown in Fig* 

6.1. Region I is the core region (U-233 plasma) . The pressure is main- 
tained at 200 atmospheres and the average core^ temperature is 25,000®K. 
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Table 6,1 Selected Properties of Moderators 


Moderator 


Atoms or 
Molecules /cm^ 

Ea(cni 

Ts(cm 

Graphite 

1.60 

0.08023x10^^^ 

2.728x10"'^ 

0.3851 

Be 

1.85 

0.1236X102** 

1.17 4x10“ 3 

0.8652 

BeO 

2.96 

0.07127x102^ 

6.771x10“*^ 

0.4846 

D 20 

1.105 

0.03323x102*^ 

3.323x10“5 

0.4519 



He (40cm) 

/Moderator 

/(20cm) 


Salt 

(bOcM) 


JIODERATOR 

(50cm) 


/ Salt Composition: 
71.7mole % LiF 
16 % BeF2 
12.3 % ThFzj 

LI COMPONENT OF SALT IS 
ENRICHED TO 99,9995% Lr , 


Fig. 6.1 Geometry Used For Moderator Comparison 





Region II is the helium region x^ith a pressure of 200 atmospheres and 
average temperature of SOOO^K* This region serves to contain the plasma 
core as well as remove heat from the plasma core. This region was main- 
tained at a thickness of 40 cm# in these initial calculations. Region III 
is the moderator region of 20 cm. thickness. This region is necessary to 
moderate the plasma core sufficiently to attain a small critical mass. 
Region IV is the molten salt region (50 cm.). The molten salt composition 
is 71.7 mole % Li F, 16 mole % BeF^, and 12.3 mole % ThF^. The Li is en- 
riched to 99.99% in the Li^ isotope in order to optimize neutron economy. 
The average temperature in this region is assumed to be 650^C. Region V 
was another moderator region (50 cm. ) . This region is necessary to main- 
tain a breeding ratio greater than 1.0. 

Results for the moderator comparison are shown in Table 6.2. It is 
evident from these calculations that either beryllium or beryllium-oxide 
could be used to meet the two primary objectives of the nuclear design. 
Beryllium- oxide was selected ma-inly due to its good thermal properties 
(m.p. of 2550^0) and also its good moderating and nuclear properties. 
Metallic beryllium could have been used if sufficient cooling were pro- 
vided to ensure no melting would occur. All subsequent calculations were 
performed using BeO as the moderator. 

Following the initial calculations to determine the most suitable 
moderating medium, calculations were performed to determine the effect 
of salt thickness, inner BeO liner thickness, and outer BeO thickness 
on critical mass and breeding ratio. These results are presented in 
Tables 6.3, 6.4 and 6.5. All dimensions except the one being studied 
remain as shown in Figure 6.1. 
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Table 6.2 Breeding Ratios and Control Masses for Various Moderators 


Moderator 


Breeding Ratio 


Critical Mass (Kg) 


Graphite 

1.06 

2900 

Be 

1.12 

14 

BeO 

1.14 

17 

D 2 O 

1.13 

80 
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Table 6.3 Effect of Salt Thickness (BeO Moderation) 


Salt Thickness 

B. R. 

Mass (KG) 

50 cm. 

1.14 

17 

40 cm. 

1.139 

16.53 

30 cm. 

1.130 

16.49 

20 cm. 

1.105 

16.33 

10 cm. 

1.026 

14.94 


Table 6*4 Effect of Inner Moderator Thickness 


Thickness B* R. Mass (KG*) 


10 cm. 

1.143 

1250 

15 cm. 

1.181 

50.5 

20 cm. 

1.140 

17 

25 cm. 

1.086 

9.5 

30 cm. 

1.024 

6.7 
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Table 6«5 Effect of Outer BeO Thickness 


Thickness 

Critical Mass (KG) 

B. R. 

50 cm. 

17 

1,14 

40 cm. 

17 

1.1423 

30 cm. 

17 

1.1420 

20 cm. 

17 

1.1414 

10 cm. 

17 

1.1397 

5 cm. 

17 

1.1376 

1 cm. 

17 

1. 13-34 
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The final reactor configuration is shown in Figure 6.2* For this 

configuration the critical core radius was 82.5 cm. The critical U-233 

mass is 26.3 Kg. and the reactor breeding ratio is 1.0099. The helium 

region is 60 cm. in thickness which gives a — ■ — = 0.727. This is 
^ ^ r-core 

sufficient for stable plasma confinement. 

The final configuration is the only one not feasible. The overall 
critical mass could be decreased but the molten salt volume would need 
to be increased to maintain a constant breeding ratio. Table 6.6 sum- 
marizes the operating parameters for the plasma core breeder reactor. 


ORIGINAL PAGE IS 
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Table 6.6 Operating Parameters for Plasma Core Breeding Reactor 


1. 

Power 

2000 MW(th) 

2. 

• kw ft 

Power Density ( /£) in Plasma 

848.8 

3. 

Core Volume (cm^) 

2.35 X 10® 

4. 

Core Radius (cm) 

82.5 

5. 

Average Thermal Plux in Plasma (^/cm^-sec) 

3.42 X 10^® 

6. 

Average Fission Density in Plasma 

2.63 X 10^3 


, fissions , q . 

C /cm^-sec) 


7. 

Critical Mass U-233 (KGS.) 

26.36 

8. 

Reactor Breeding Ratio 

1.0099 

9. 

Peak/ Average Fission Density in Plasma 

1.126 

10. 

Fuel Absorption/Fissions in Plasma 

1.113 
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7. DESIGN CONSIDERATIONS OF THE BREEDER REACTOR 


The material chosen for the pressure vessel was type 347 stainless 
steel. The pressure vessel is spherical and calculated minimum pressure 
vessel thicknesses are indicated" on Figure 7.1. The designed operating 
temperature of the pressure vessel is 422°K. The inside radius is 193 cm. 
with an operating pressure of 200 atm. The minimum wall thickness is 
then 15.3 cm. and the designed wall thickness is 20 cm. 

Figure 7 . 2 illustrates the helium and molten salt flow paths for the 
breeder reactor. Approximately 7% of the fission energy appears as energy 
of the neutrons and gamma rays. This energy, 140 MW, will ultimately be 
deposited in the reactor moderator, molten salt, and other structural mate- 
rials. Also, 0.1% of the radiated energy from the plasma core; or 1.86 MW, 
will be deposited at the BeO wall- of the cavity. This heat will be absorbed 
by the helium passing through the porous wall. The helium cooling the molten 
salt enters the heat exchanger at 422°K and exits at 649°K. 

The molten salt enters the reactor at 920*^K (nearly the minimum tem- 
perature possible) and is directed by Zircaloy baffles radially inward. 

The molten salt .cools the inner BeO moderator layer and then leaves the 
reactor at lllO^K. Assuming that the salt absorbs all 140 MW of the heat 
from the inner layer, the flow rate of the salt will be 542 kg/sec and 
assuming a 10% pressure loss in the molten salt system, the needed pump- 
ing work is 330 KW. Part of the molten salt from the hot leg, 0.332 gm/sec, 
is diverted to the reprocessing systems. The amount of salt in the 
reactor is 1.75x10^ kg and the time spent by the salt in the reactor is 
32 sec. There is a 2 cm. gap filled with helium between the molten salt 
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and the outer BeO moderator layer. This is to insulate the outer modera- 
tor and the pressure vessel from the high temperature molten salt. 

Zircaloy is used for internal structural materials because of its 
low neutron absorption. There is a thin Zircaloy layer between the 
various regions of the reactor. The inner porous liner is a Zircaloy 
shell filled with BeO, 25% of the volume is passages for the helium coolant. 
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8. REPROCESSING SYSTEM FOR PLASMA CORE BREEDER REACTOR 


Since it is desirable to have the Plasma Core Breeder Reactor (PCBR) 
be a self “-contained unit, generating its o\m new fuel, an on-^line repro- 
cessing system for the molten salt blanket is a necessity. This chapter 
describes protactinium removal and salt purification processes, calcula- 
tions of expected flow rates, and equilibrium concentrations of various 
isotopes present in the system, 

The salt used in the blanket is an eutectic mixture composed of LiF, 
BeF 2 , and ThF^ in the ratios of 72:16:12 mole percent. This particular 
combination was developed at Oak Ridge National Lab in conjunction with 
the Molten Salt Breeder Reactor program. 

When thorium atoms contained in the salt enter the neutron field 
around the core, some of the atcms absorb a neutron and transmute to pro- 
tactinium as shown in Figure 8.1. The protactinium eventually decays to 
uranium which can then be fed to the core as new fuel. However, as seen 
in Figure 8.1, Pa^^^ has a substantial cross section (22 barns) and since 
its half life is 27 days, Pa acts as a poison, siphoning off neutrons 
which could otherwise irradiate Th atoms. For this reason, it is desir- 
able to remove Pa from the molten salt loop and allow it to decay outside 
the core. 

However, since it is impossible to have a zero protactinium concen- 
tration in the molten salt blanket, there will be some uranium present 
in the core. Some of these atoms will fission and, consequently, there 
will be some uranium fission products in the molten salt loop. Some of 
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Fig. 8.1 The chain of isotopes created by neutron irradiation 
of Th232. 
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these fission products have large cross sections as shown in Table 8.1. 
Note that Xe and other gaseous fission product poisons are not listed be- 
cause it is assumed that the blanket can be vented and these gaseous 
products easily removed. As will be shown later, the necessity of keep- ' 
ing the concentration of fission products at a low level determines the 
amount of time which the salt can stay in the irradiated blanket region. 

In order to achieve the above neutronics goals, a fluorinat ion-reduc- 
tive extraction system was developed at Oak Ridge National Lab. A des- 

2 

cription of this process is as follows: 

The fluorinat ion-reductive extraction system for isolating 
protactinium is shown in its simplest form in Figure 8.2. 

The salt stream from the reactor first passes through a 
fluorinator, where most of the uranium is removed by fluor- 
ination. Approximately 90% of the salt leaving the fluor- 
inator is fed to an extraction column,* where it is counter- 
currently contacted with a bismuth stream containing lithium 
and thorium. The uranium is preferentially removed from 
the salt in the lower extractor, and the protactinium is re- 
moved by the upper contactor. A tank through which the bismuth 
flows is provided for retaining most of ^ the protactinium in 
the system. 

The bismuth stream leaving the lower contactor contains 
some protactinium as well as the uranium that was not removed 
in the fluorinator and the uranium that was produced by the 
decay of protactinium. This stream is contacted with a 
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Table 8.1 Rare Earth Fission Product Absorption Cross Section 




Nd-143 

330 barns 

La-139 

8*9 barns 

Eu-153 

320 barns 
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Fig. 8.2 Plasma Core Reactor Reprocessing System 
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mixture in the presence of approximately 10 % of the salt 
leaving the flourinator in order to transfer the uranium 
and the protactinium to the salt. The salt stream, contain-- 
ing UF 14 and PaFi,., is then returned to a point upstream of 
the fluorinator, where most of the uranium is removed. The 
protactinium passes through the fluorinator and is subse- 
quently extracted into the bismuth. Reduct ant (Li and Th) 
is added to the Bi stream leaving the oxidizer, and the re- 
sulting stream is returned to the upper contractor. The 
salt stream leaving the upper contactor is essentially free 
of uranium and protactinium and would be processed (for 
removal of any fission product gases and additional thorium 
added to compensate for that which had been consumed.) 

Figure 8.3 describes the UF. to U metal conversion process, Unfor- 

0 

tunately this is a batch process instead of a continuous flow system 
as is present in the remainder of the reprocessing set-up. However, 
there should be no problem providing temporary storage tanks for UF^. 

The IIF 5 initially enters a reaction chamber where it is mixed with 
hydrogen. A reaction is triggered and UF 4 powder and HF gas is produced. 
The UF 4 is then loaded into a steel '^bomb” which has been coated x^ith 
fused dolomitic lime — lime is one of the few oxides that does not react 
with molten uranium. The ”bomb” is then heated to 565® C where an exo- 
thermic reaction takes place and uranium metal solidifies on the bottom 

of the "bomb”. The MgF£ is removed and U metal of high purity can then 

4 

be taken from the bottom of the "bomb" and sent to the reactor. 


C'> 
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Process 



Given certain constraints on the reprocessing system it is possi- 
ble to calculate the flow rates which would exist in both the molten 
salt and bismuth loops* It is also possible to calculate protactinium 
concentrations throughout the reprocessing system and therefore deter- 
mine uranium concentrations throughout the system. The constraints 
which are placed on the reprocessing system are as follows: 

1) The protactinitim concentration in the molten salt blanket is 
allowed to reach 95% of the equilibrium value obtained if the salt re- 
mained in the active region of the reactor for an infinite amount of 
time, provided that the concentration of protactinium does not cause 
parasitic absorption of neutrons by fission products greater than 1% of 
the absorptions which are due to thorium captures. 

2) The volume of the blanket and the flux in the blanket shall be 
determined by breeding ratio constraints as explained elsewhere in this 
report. 

3) The uranium removal efficiency of the f luorinator and oxidizer 
is 98%.^ 

4) The operating temperature of the system shall be 640°C (neces- 

4 w 

sary because the salt is a eutectic mixture) . 

5) The Li concentration in the Bi loop shall be 1%. The Th con- 
centration in the Bi loop shall be held at less than 50% of the solubil- 

5 

ity of Th in Bi. 

6) The Pa distribution coefficient for the contactors, defined as 
(mole fraction of Pa in Bi at equilibrium) / (mole fraction of Pa in salt 

5 

at equilibrium), is taken to be 100. 
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7) The following physics data is assumed: 


Neutron Flux 
Volume of Blanket 
Molar Volume of Salt 
Molar Volume of Bi 
Pa Absorption Cross section 
Th Absorption Cross section 
U Absorption Cross section 
U Fission Cross section 
Ea Decay Constant 


l.OE+14 n/cc~sec. 
5.29E+06 cc. 
0,0598 moles/ cc.. 
0.0469 moies/cc. 
22 barns 
7.4 barns 
571.3 barns 
524.5 barns 
2.97E-07 sec“l 
4.32E+21 


Concentration of Th in salt 
8) Due to its very short half-life Th^^^ is assumed to transmute 
directly to upon being struck with a neutron* 

To satisfy assumption 1, we must examine if the Pa concentration in 
the salt from the output of the blanket will be governed by the rate of 


fission product captures. To determine the number of fission product 
captures we must first solve for the Pa and U concentrations. This is 
done as follows: 


+ APa = a ^ Th (8.1) 

dt ^ 

where is the flux, Th is the thorium concentration, and A the Pa decay 
constant. 

Solving Eq. 8.1 gives 


Pa 



-At 

e 


o ^^<f) Th 
a 

_ 


-) 


Pa 

o 


( 8 . 2 ) 
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The equation for the uranium concentration as a function of time is 


“ = - (j) 0 “ U + X Pa (8.3) 

nt a 

where U is the U-233 concentration. 

Solving this equation we have 



If a given atom of material is assumed to spend time T in the blan^ 
ket, then the number of fissions which occurs during this time is 


No. of fissions 


= ] 0^Vu(t)dt 


(8.5) 


Evaluating this integral we have 


u 


No. of fissions = 


Th , 

Th , 

a ® 


-0 

a ^ 


u 




(0 “<}.-X) 


m ,m 

Th^ 

r , -0 4 >T 


0, <^Th 

, -XT a 


a p 

1-e e - 1 


X ao 

L ^ 




- 


( 8 . 6 ) 


and the fission product concentration at the end of a cycle of length T 
is given by 


[F.P.] = 


pT (f,t 

\ a4> U(t) e dt 

‘'o 


U 

— 0£ ({)T 

e «y (Wo. of fissions) 


(8.7) 
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where y is the probability per fission of getting a particular fission 
product* Since the f luonixater removes 98% of the uranium in the molten 
salt on each pass through the system, the entering concentration to the 
blanket region can be taken as effectively zero. 

Solving Eq. 8,7 for a variety of times T, the results, given 

^Eu 

where is the absorption cross section of one of the most 

troublesome rare earth fission products, Eu^^^, are shoxm in Table 8.2. 
It should be stated that this estimate of the Eu^^^ concentration is 
high due to the approximation in Eq. 8.7. However, even with this 
high estimate it can be seen that no fission product removal system is 
necessary. 

To determine the flow rates and concentrations in the system, one 
must make use of the following mass balance equations.^ If one refers 
to the hypothetical exchange column shown in Figure 8.4 



Figure 8.4: Exchange Column Flows 

then a material balance yields the following equation: 


or 

^ " V (y^ y^) (8.9) 
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Table 8.2 Protactiniiun and Europivim Concentrations in Blanket 


Time in Blanket Pa Concentration 

(days) (% of equilibrium) 



114 

.95 

6 

45.5 

.70 

.75 

26.1 

.50 

.165 

19.1 

.40 

.072 

10.7 

.25 

.0135 

3.91 

.10 

6.9(10' 
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where L and V are flow rates in moles /sec and x and y are concentrations 
of the transferring material expressed in mole fractions. Now at equili- ' 
britim 

Yl K - x^ (8.10) 

where K is a constant known as the distribution coefficient. Substitu*- 
ting for x^ in Equation 8.9 and solving for we have 


y 


1 


= 




X 

o 


( 8 . 11 ) 


So if we knew the two inlet concentrations and if we can find the flow 
rates then the outlet concentrations can be calculated. 

The value of the flow rates in the Bi and blanket loops must be 
solved for iteratively. A flow chart of the solution process is shown 
in Fig. 8.5 and a FORTRAN program written to solve this problem. A 
value for the Bi flow rate is assumed and for given Pa core concentra- 
tion, neutron flux, and core volume, the flow rate in the blanket, res- 
idence time in the core, and input concentration of Pa to the core can 
be solved for iteratively. 

Now, as referenced above, ORNL report number 4344 gives the distri- 
bution coefficient of Pa as a function of time of contact and relative 


volumes of salt and Bi. Picking a specific distribution coefficient 
determines the time of contact and the relative volume of the two com- 
ponents. A new value for the Bi flow rate can then be calculated by 

A 

using the value of the blanket flow rate calculated above. The entire 
iterative procedure is then repeated with the new Bi flow rate. 
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Set Pe 2» Flux 

i 


Assume Pai, Pa3, Bi flow rate 

■ I 

-^Calculate- Time = fn(Pa2» Pai, Flux) 

i ■ 

Calculate Core flow rate = fn(Time) 


i 

-Calculate Pai = fn(Pa2, Core flow, Bi flow, Pas) 

(when converge) 


-Calculate Bi flow rate = fn(all variables) 

(when converge) 
Calculate Pa4 = fn(Pai) 

i 

Calculate Pa3 = fnCPati) 

(when converge) 

Stop 


Pai = Core input Pa concentration 

Pa2 = Core output Pa concentration 

Paa = Bi loop contactor input Pa concentration 

Pa4 = Bi loop contactor output Pa concentration 


Figure 8.5 Flowchart for Calculation of Reprocessing 
System Flow Rates and Pa Concentration 
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Once the flow rates have been calculated, the output Fa concentra-- 
tion in the Bi loop from the contactor can then be found from Eq. 8.11 
and the input concentration from Eq, 8.9. 

It should be noted at this point that if a contactor is composed of 
several stages with K being the distribution coefficient in each stage, 
then the procedure described above can be applied to the whole system 
with the number of. stages, N, given by the expression^ 


N = 


log 


rH' 

1 

f^n+l ■ 


[A 1 

i 

^ A 


( 8 . 12 ) 


log A 

where A is the absorption factor and is defined by A = L/(KV), However, 
for simplicity, this report assumes the contactor to have only one stage. 

Performing the calculations described above, the reprocessing system 
parameters shown in Table 8.3 were calculated. 

In conclusion, an on-line, on-site reprocessing system for the Plasma 
Core Reactor is quite feasible. The technology is available today and the 
chemical processes involved in uranium separation have been proven in vari- 
ous laboratory experiments at Oak Ridge National Lab as a part of the Molten 
Salt Breeder program. Only a Pa isolation system is required, no fission 
product removal system is necessary. 
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Table 8.3 


Summary of Reprocessing Systems Parameters for /PCBR 

/ 

The fission concentration in the blanket per cycle is .333E+18 fissions 
per cc. 

The effective time spent in the core is .655E+01 days. 

The flow rate through the core is .691E+00 G-moles/sec. 

Input protactinium concentration to the core is .549E+18 atoms per cc. 
Output protactinium concentration from the core is .145E+20 atoms per cc. 
The flow rate in the Bi loop is .180E+00 G~mole/sec. 

The Pa concentration in the Bi loop entering the contactor is .848E+18 
atoms per cc. 

Output Pa concentration from the contactor is .430E+20 atoms per cc. 

The number of stages in the contactor “is“.100E+01. 
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9. MHD GENERATOR AND SEPARATOR 


The plasma core reactor-MHD system was first explored by Colgate 
and Aamodt in 1957.^ In 1973, Williams and Clement presented calculations- 
for plasma core reactor-MHD power plants which had efficiencies of 70 
percent or more." 

The problems associated with the plasma core reactor-MHD system are 

likely to be severe, but the potential of this system is large enough to 

warrant serious investigation. The growing program in MHD power genera- 

3 

tion and the ongoing UF. reactor tests at Los Alamos will provide 

o 

information related to some critical components. However, studies are 

needed to define the problems unique to plasma core reactor-MHD systems, 

to offer possible solutions, and to formulate an experimental program 

if such a program is desired. In view of the growing interests in MHD 

and plasma core reactor technologies and long lead times in research 

and development of both space and terrestrial applications, it is both 

timely and beneficial to initiate such studies. 

The analysis of the MHD generator follows that of Ref. 2. The 

generator is a segmented electrode Faraday generator with cesium seeded 

helium as the working fluid. Since He does not ionize significantly 

until SOOO^K, even at 1.013 pascals (10^^ atmospheres), it may be treated 

4 5 

as a perfect gas. * The relevant gas properties are listed in Table 9.1. 


Table 9.1 Gas Properties of Helium 


Ratio of Specific Heats, y 

Specific Heat at Constant Pressure, Cp 

Gas Constant 


1. 6667 

5192 ~— 
Kg°K 

2077 

Kg®K 
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The electrical conductivity of the seeded gas is a function of 
temperature, pressure, and seed mass fraction. Only thermal ionization 
was considered in this study, although the actual generator includes 
several non- equilibrium processes due to slow recombination during 
expansion through the duct (frozen flow) , current caused by a motional 
emf or an electric field, fission fragment, and electromagnetic radiation. 
The design parameter that is affected directly by the electrical con- 
ductivity is the length of the generator. This in turn will influence 
other parameters as will be subsequently seen. 

The generator is a constant velocity generator which is divided 
into 15 segments to eliminate the Hall current. The state points which 
are used in the analysis are shown on Fig. 9.1. Note that for a unit mass 
of helium, a fraction x is diverted to cool the nozzle, a fraction y is 
diverted to cool the electrodes, and the remaining fraction 1-x-y cools 
the blanket before entering the reactor cavity. 

Given the cavity power, Q^, and the cavity inlet and exit stagnation 

temperatures, T and T , respectively, the mass flow rate at the cavity 
tl t£ 

exit is 


= H — 

t2 tl 




(9.1) 


where H aiuJ H are stagnation enthalpies at the reactor entrance and 
t2 tl 

exit, respectively. 

The static temperature and pressure at the nozzle exit (MHD duct 
entrance) for an isentropic process are 


T 


Ts 


t2 


Y~ 1 2 

1 + ^ Ms 


(9.2) 


94 


OB' POOR 



1 kg/ sec 



"Plasma core 
reactor 


Nozzle 



Uranium 

processing 


Y-1 2 

1 + -V" Mg Y-1 


(9.3) 


where Mg, the Mach number at the MHD entrance and stagnation 

pressure at the cavity exit are input quantities. 

The helium velocity at the MHD inlet is 


Dg = V yR Tg Mg 


This velocity is constant across the generator. The kinetic energy 
of the gas is 


K.E. = 


The total enthalpy at the nozzle exit has a component from the cavity 

and a component due to nozzle cooling 

(l-x-y)H + X H 

*2 t. 


The static enthalpy at the MHD inlet is 


Hg = H - K.E. 

tg 


The corresponding static pressure and temperature at the MHD inlet are 
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P3 

= P3 

(9.8) 

T 3 

H 3 

(9.9) 


The inlet density of the gas at the MHD inlet is 


P3 - 


P3 


(9.10) 


while the mass flow rate is 


» • 


1-y 

m = mo - — — 
^ 1-x-y 


(9.11) 


The MHD inlet area is then 


A3 “ 


m3 


P3U3 


( 9 . 12 ) 


The pressure ratio for each segment is 

1 





n 


P3 


(9.13) 


where n is the number of segments and , the pressure ratio 

Vk 

across the generator are given. 

The exit pressure at each segment is 


'm 


i+1 


'm. 

1 


m 


X“ij. . n 


(9.14) 


with p^ = P3. The pressure drop for each segment is 


Ap. = p - p 

x-hl X 


X ' 1 ^ ^ 


(9.15) 
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th 

The exit temperature of the i segment neglecting electrode 


cooling is 


T 

“i +1 


T 

m. 


m 


1 

K(y-1) 


Y 


1 — 1^* • * y n 


(9.16) 


where K is the MHD loading factor and T ^ = Ta. 

mi 

The exit enthalpy of each segment with transpirational cooling 
of the electrode is 


H 

m 


H 


i +1 


[ 1 - y + (i-1) -^1 + H 


^i+1 


(1 - y + 1 -^) 
n 


i = 1, . * . yXi 


(9.17) 


with H = = H 3 . The exit temperature at each segment is then 


mi 


H 

m 


m 


i+1 


i +1 

Cp 


i 1^9 m « ^ XI. 


(9,18) 


with “Is- lo find the average conductivity in each segment, cr^, 
requires the average temperature and pressure in each segment given by 


T 


m^ 

1 


T + T 

”i+l 

2 


(9.19) 





2 


(9.20) 


The length of each segment is given by 
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AL. = 

X 


B^TJ3d^(l-K) 


where B is the magnetic flux density and is an input quantity. 

The density corresponding to T and p is. 

i X 


m. 

1 


_n^ 

RT 

m. 


X'^1^. « .^n 


which is used to determine the inlet area of each segment 

*3 (i-1) 

mo + 

A „ n 1-v 

X = 1,. . . ,n 


Us 


m. 

- X 


The exit area of each segment is 




X P 

m Y -1 


i = 1, . - . ,n 


The generator length is 


n 

L = r AL. 

««*• -I 

x=i ^ 


“ Pjjj and so that the MHD exit enthalpy is 


n+1 


n+1 


= Cp T^ 


(9.21) 


(9.22) 


(9,23) 


(9.24) 


(9.25) 


(9.26) 
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The stagnation enthalpy at the MHD exit is 


H 


t 


4 


H. 



The MHD exit Mach number is 

U3 

M = 

4 YyO!;- 


(9.27) 


(9.28) 


After leaving the MHD duct, the Mach number of helium is reduced to 
0.1 by a diffuser before entering the separators. The temperature and 
pressure at the separator exit are 


T, = T, 
5 4 


(H- ) 

^ «5^> 


(9.29) 


iY-2- 


P5 = P4 


r (1 + ^ ) < 11 

[ (1 + ^ J 


(9.30) 


Assuming no losses, the thermal energy in the MHD generator is 
equal to the electric power produced. 


QiAur. = m^(H^. +~ — H^ + -r-^ — H^. -- — ~ H. ) 
^MHD 2 '- t^ 1-x-y 1-x-y 1-x-y ' 


(9.31) 


Since two separators are used, the flow rate into each separator is 
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s, 

in 


(9,32) 


in 


1 “3 


l-y 2 


The gas velocity at the separator exit is 


Ug = VyR Tg Mg 


(9.33) 


The helium density at the inlet or exit of each separator is 

P5 

^5 ~ R Tc 


(9.34) 


The inlet area of each separator is 

A 


m. 

in 


(9.35) 


in P 5 U 5 

The separator has two exits connected with two turbine-compressor 


units. The exit area is 


m 


ex 


ex 


P U 
^5 5 


(9.36) 


where 


m 


ex 


m. 


xn 


(9.37) 


A computer program was developed for the MHD generator and sepa- 
rator. A baseline case was established and is summarized in Table 9.2. 
The assumed magnetic flux density of 18.0 Teslas is very high compared 
to fields of normal MHD generators but is within superconducting 
magnet technology. This large value was due to four constraints on 
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Table 9.2 MHD Reference Design 


Core Power = 2000 JlWt 

MHD Electric Power = 1022 MWe 

Reactor Exit Temperature = 3500°K 

Cavity Pressure = 2.027x10^ Pascals (200 ATM.) 

MHD Pressure Ratio, = 3.0 

MHD Exit Temperature = 2273”K 

MHD Inlet Mach Number = 0.500 

MHD Exit Mach Number = 0.596 

Helium Gas Velocity in MHD Duct = 1672 M/sec. 

Load Factor = 0.8 

Magnetic Flux Density = 18.0 Teslas 

MHD Length = 3.47 M 

Length/Diameter = 10.75 

MHD Inlet Diameter = 0.222 M 

MHD Exit Diameter = 0.323 M 

CS Seed Fraction = 0.015 

Number of Electrode Segments = 15 

Conductivity in First Segment = 206 MHOS/M 

Conductivity in Last Segment = 7.1 MHOS/M 
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the generator. First, a reasonable length had to be found. A large 
pressure ratio would result in excessive lengths (greater than 10 meters) . 
The length of the generator was a strong function of the pressure ratio 
(Fig. 9.2). This effect was due to the strong variation of electrical 
conductivity with temperature. At low exit temperatures (less than 
2000^K), corresponding to a large pressure ratio, the conductivity de- 
creases by several orders of magnitude from its value at the MHD inlet 
resulting in a very long generator. 

Even more important was the length to diameter ratio of the generator 
as a function of the MHD pressure ratio (Fig. 9.2). For inert gas genera- 
tors, L/D is about 10 to insure a well behaved boundary layer. This 
value may be conservative but was taken as a constraint on the system. 
Figure 9.3 shows that the magnetic flux density is also a strong function 
of the MHD pressure ratio for a L/D of 10. 

However, the MHD pressure ratio cannot be taken to be too low; other- 
wise, very little power is extracted. It is desirable to convert as 
much of the thermal energy in the fluid in the MHD duct as this leads 
to a higher plant efficiency. It was desired to extract 1000 MWe or 
more from the plant at high plant efficiencies. This factor plus the 
constraints on length and L/D led to the values in Table 9.2, 

The magnetic field may be considered too high. In this case, the 
MHD pressure ratio can be dropped which lowers the magnetic field 
requirements (Fig. 9.3), but decreases the electrical power output (Fig. 
9.4). Another alternative is to keep the electric output constant, but 
increase the core power and decrease the pressure ratio. Again the 
magnetic field requirement is lower, but the efficiency of the plant 
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1 ' 


32 



2 *■ 

^1*5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

MHD Pressure Ratio 


Fig. 9.3 Magnetic Flux Density Versus MHD Pressure 
•Ratio for L/D = 10 
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decreases* The selection of a proper power level depends on economic 
factors such as capital and fuel costs which is beyond the scope of 
this study* 

In addition, operating at lower cavity pressures would result in 
an increase in electrical conductivity as large pressures inhibits 
ionization* However, 200 atmospheres was chosen as the cavity pressure 
to insure low critical masses for the breeder and actinide transmutation 
reactor concepts* No optimization of pressure was performed for the 
combined reactor-MHD system* 

The constraints of electric power output, length, length-^to-niiame- 
ter ratio, and cavity pressure lead to a MHD generator with a large 
magnetic field but is within the technology of ^ superconducting magnets* 

An optimization study is needed on cavity power and cavity pressure 
to obtain a better MHD generator* Also, a two temperature model is needed 
for the calculation of the electrical conductivity to account for non- 
equilibrium processes* 

The main question of the MHD generator that will have to be answered 
or solved concerns the flow of uranitim through the duct* The fission 
fragments may enhance ionization in the channel but may also cause 
serious problems. Some of the problems are: 

a) Shielding the superconducting magnet from neutrons and gammas* 

This should not be too much of a technical problem but may be an 
economic one* Research on fusion reactors should provide some 
information in this area, 

b) Radiation damage to the electrodes over a prolong period of time. 
Not much is known in this area and it should be given attention, 
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c) Condensation of uranium droplets may short out- the electrodes. 
This is probably the most important concern of this type of 
system and it deserves considerable research. 
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10. POWER PLANT SYSTEMS 


In order to achieve maximum effectiveness from the high temperature 
coolants from plasma core reactors ^ it was decided that a ternary 
power cycle would produce the highest efficiency power plant. The ter- 
nary cycle consists of a combination of MHD, gas turbine, and Rankine 
cycle energy conversion units. Two concepts x^ere considered — systems 
with and without a high temperature regenerator in the helium loop. 

The objectives of the study were as follows: 

1. Model the nuclear MHD power plant cycle. 

2. Analyze the power output from the three energy conversion units 

and evaluate plant overall efficiency. 

3. Make a parametric study of the effect of changing operating 

variables on plant overall performance. 

The components of the propsed systems are shown in Figs. 10.1 and 
10.2. For both systems the heat source is a high temperature reactor 
with a uranium plasma core (UPC) . Helium is used to cool the core and is 
the working fluid for the MHD generator and gas turbine. Helium enters 
the MHD generator and expands from state 1 to state 2. After the MHD 
channel it passes through a diffuser and enters the separator where the 
uranium is separated from the helium. Because of its high temperature 
at the exit of the separator, the helium needs to be cooled. This process 
takes place in a mixing tank (MT) where high temperature helium is 
mixed with helium coming from compressor (Cl) in Fig. 10.1. At point 4 
the cooled helium enters the gas turbine. After expansion, helium passes 
through heat exchangers (HEl) and (HE2) where heat is transferred to the 
steam cycle. Helium is then compressed in the first stage (Cl) and 
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Fig. 10.2 Nuclear MED Power Plant With Regeneration 



Plasma Core; MHD *- MHD Generator; SEP Diffuser and Separator; HEi — Regenerative 
Heat Exchanger; HE2,3,4 and 5 - Heat Exchangers; GE - Electric Generator; GT - Gas Turbine; 

Compressor; ME - Electric Motor; ST - Steam Turbine; COND - Condenser; WP - Water Pump; 

FWH - Peed Water Heater, 









then split into two loops. The inner loop is provided to feed the 
mixing tank. The outer loop is to cool the reactor and for this purpose 
the helium needs two more compression stages (C2 and C3) . The inter- 
coolers (HE3 and HE4) transfers heat to the steam cycle. Before enter- 
ing the reactor, two fractions, x and y, are taken for cooling the nozzle 
and MHD duct. Since the MHD pressure ratio is greater than that of the 
gas turbine, and to maintain a functional unity for the components with 
the same mass flow rate of gas, compressor (Cl) is on the same axis with 
the gas turbine and electric generator (GE2) , while compressors (C2) and 
(C3) are powered by electric motors. 

The nuclear MHD power plant with regeneration (Fig. 10.2) is similar 
to the previous system in the heat source region (reactor, MHD, and SEP). 
After the separator, helium passes through a regenerative heat enchanger 
(HEl) and enters the gas turbine. Between states 5 (exit GT) and 12 
(inlet breeding zone of the reactor), helium is compressed in three 
stages and transfers a part of its heat to the steam cycle in heat ex- 
changers HE2, HE3, HE4, and HE5. Unlike the previous system, helium passes 
through only one loop and transfers a major part of its heat through 
regeneration. 

For both systems the steam cycle is the same. Water heated in one 
or two stages of feedwater heaters (or directly from the condenser) goes 
to a boiler. Superheated steam enters a steam turbine (ST) and then • 
passes through a condenser (COND) . 

All studies used values for input data according to current com- 
mercial technology (i.e. efficiencies for steam cycle components, gas 
turbine, and compressors) or with current use in MHD research. For the 
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ORIGINAL RAGS IS 
OF POOR QUALITY 



isentropic efficiency of the MHD generator we assumed values of 75% - 
1 2 

80%. ’ The electric efficiency of the MHD generator was provided by 
Fig. 10.3. 

The modeling of the MHD cycle consisted of defining a pseudo - 
Brayton cycle and treating the expansion within the MHD generator in 
a similar manner as in a gas turbine. In order to analyze fhe two 
systems it was necessary to write two computer codes: 

(1) HMHD-1 — code to analyze the nuclear MHD power plant without 

regeneration in the helium loop 

(2) NMRD-2 — code to analyze the nuclear MHD power plant with 

regeneration in the helium loop. 

The basic logic followed in the computer codes is presented in 
Fig. 10,4. Table 10.1 lists input parameters for each system. 

The codes are general in that they permit any changes in input data. 
From the input data, using special subroutines (TSAT, SUPER and SATL) , 
the program STEAM evaluates all necessary parameters for the steam cycle 
and calculates the net power produced within the cycle. In addition, 
the code prepares the enthalpy values for determining the mass flow rate 
ratio between the helium cycle and steam cycle. For evaluation of state 
parameters characteristic to the top cycles the codes have implemented 
a subprogram MHD (different for the two codes) . Taking information 
from calculations done by STEAM and MHD the codes evaluate the power 
distribution for each energy conversion unit and calculates the plant 
overall efficiencies. 

For a pair of selected parameters by the user, the codes permit 
a parametric study of the whole system yielding information for evalu- 
ating power distributions and overall efficiencies. The first parameter 
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Call Subprogram Steam (Boiler Exit 
Temperature and Pressure) . Eval- 
uate steam cycle parameters and 
return: steam work/lb steam, 
boiler entrance temperature and 
enthalpy. 


Call Subprogram MHD (input data, 
boiler entrance temperature and 
enthalpy) . Evaluate He-cycle 
parameters and return work MHD, 
work compressors, He-steam ratio. 



Calculate power distribution and 
plant overall efficiency. 


* 

Increment first ; 
ter (20 times) 

Indicated parame- 



Increment second 
ter (4 times) 

indicated parame- 


Print out result 
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Table 10.1 Input Data for NMHD-1 and NMHD-2 


Index 

NMHD-1 

NMHD-2 

1 

Boiler Temperature 

1000°F 

Boiler Temperature 

1000. ®F 

2 

Boiler Pressure 

1600 psia 

Boiler Pressure 

1600 psia 

3 

Condenser Pressure 

1*0 psia 

Condenser Pressure 

1.0 psia 

4 

Steam Turbine Efficiency 

81% 

Steam Turbine Efficiency 

81% 

5 

Pump Efficiency 

80% 

Number of Feed Heaters 

0,1 or 2 

6 

Number of Feed Heaters 

0,1 or 2 

Reactor Temp Difference 

200°K 

7 

Compressor Efficiency - 

85% 

Compressor Efficiency - 

85% 

8 

MHD Inlet Temp 

3000°K 

MHD Inlet Temp 

3000°K 

9 

MHD Inlet Press 

200 bar 

MHD Inlet Press ■' 

200 bar 

10 

MHD Pressure Ratio 

5.0 

MHD Pressure Ratio 

3.0 

11 

Gas Turbine Pressure Ratio 

2,0 

Gas Turbine Press, ratio 

3.0 

12 

Feed Heater 1 Pressure 

12. psia 

Feed Heater 1 press. — 

12. psia 

13 

Feed Heater 2 Pressure - 

4. psia 

Feed Heater 2 press. — 

4.0 psia 

14 

Bottom Temp Difference - 

150°K 

Bottom Temp Diff * 

150°K 

15 

MHD Inlet Mach No. 

0.5 

MHD Inlet Mach No. 

0.5 

16 

Sep Outlet Mach No. 

0 ._ 1 _ 

Sep Outlet Mach No. 

0.1 

17 

Gas Turbine Inlet Temp - 

1500*^K 

Gas Turbine Inlet Temp 

1500°K 

18 

MHD Efficiency 

49% 

MHD Efficiency 

49% 

19 

Gas Turbine Efficiency - 

85% 

Gas Turbine Efficiency 

85% 

20 

Number of Compress Stages 

3.0 

Number of Compress Stages 3.0 
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is incremented twenty times; for each step of this variation, the 
second parameter is incremented four times • 

The objectives of the parametric studies were to establish the 
influence of different parameters on overall efficiencies for each system 
and to determine power distribution- For a suitable comparison of the 
influence of each parameter on the overall efficiency, a new parameter 
called ^^sensitivity” was defined - 

The sensitivity of the plant overall efficiency is defined as the 

ratio 

n- 

sensitivity (%) ” ^ (10-1) 

where: r\ is the reference value efficiency obtained in the first step 

o 

of calculation; 

th 

is the efficiency obtained in the i step of the parametric 
study - 

The results obtained in this study provide valuable information 
concerning the behaviour of overall efficiencies. Results are presented 
in Fig- 5- The most important parameters affecting overall efficiency 
are as follows: 

1- WHD inlet temperature- 
2. MHD pressure ratio - 

3- Temperature drop across the regenerative heat exchanger 
(for M1HD2). 

4. Gas turbine pressure ratio 

5- MHD inlet pressure - 

6- Boiler temperature 

7- Boiler pressure 
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Sensitivity of the Plant Overall Efficiency 



0 2 4 6 8 10 12 14 16 18 20 


Number of Increment Steps (i) 

1 MHD Inlet Temperature {3000-5000°K, =62.75%) 

2 MHD Pressure Ratio (1.5-11.5,7^=59.65%) 

3 Temp Drop Regenerator (300-1000®K, 7,=60.13%) 

4 Gas Turbine Press. Ratio (1 . 5-11 . 5 , 7o=61.45%) 

5 MHD Inlet Pressure (100-200bar , 7, =58 . 64%) 

6 Boiler Temperature ^000-2000®F, 7o=62 .75%) 

7 Boiler Pressure (1000-2000psia, 70=62. 42%) 

8 Nxmber of Comp. Stages (1 . 0-10 .0 , f*=67.5]%) 

9 Bottom Temp. Difference (100-300®F, 7^=65.10%) 


Fig. 10.5 Sensitivity of the Plant Net Overall Efficiency to 
the Variation of the Main Parameters 
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8* Number of compression stages 

9. Bottom temperature difference* 

Results from the sensitivity study are shown in Fig* 10-5. 

After performing the sensitivity analysis a study was made of the 
effect on overall efficiency of varying the reactor coolant outlet 
temperature from 3000^K to 4000^K for the two systems. Tables 10.2 
and 10.3 list the output works for each energy conversion device and 
plant overall efficiencies for systems without and with regeneration on 
the helium loop, respectively. The effect of reactor outlet temperature 
is more substantial on the system with regeneration. The effect of 
reactor outlet temperature on plant overall efficiency is shown graphi- 
cally in Fig. 10.6. 

For Nuclear MHD Power Plants without regeneration, the major contri- 
bution of electric power is due to the steam turbine subsystem (36-03% - 
36.36% from 100 % heat produced by the reactor). Due to a significant 
fraction of the electric power being produced by the steam cycle with a 
low efficiency (40%), it is desirable to shift the power production 
toward the top of the cycle to improve the overall efficiency. This can 
be achieved by reducing the mass flow rate of helium within the inner 
loop and increasing the pressure ratio of the MHD generator- This 
system produced overall efficiencies that are 15-20% higher than actual 
power plants in use and that are 5-10% higher than the expected coal"' 
fired MHD power plant. Due to the relatively low temperatures within 
the helium loop, this type of power plant could be considered as a first 
step in a national program of implementation of MHD power plants with 
a nuclear heat source. 
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Table 1Q>2 Plant. Het Overall Efficiencies For MHD Inlet Temperature Variation 


MHD Inlet 
Temperature 

j 3000 '’k 

- 

3250°K . 

3500°K 

1 3750 


4000®K 

Gas Flow Rate 
Through the GT* 

mgHH 

2*60 kg/sec 

2.88 kg/sec 

3.15 kg/sec 

3.42 kg/sec 


12265*71 

100*0% 


100.0% 

14862.21 

100*0% 

16160.46 

100.0% 



u 

' MHD 

1777.71 

14.49% 

2077*87 

15*32% 

2378.55 

16.0% 

2679.22 

16.58% 

2929.90 

17.07% 

'^GT 

456.46 

3*72% 

510*00 

3.76% 

563.54 

3.79% 

617.68 

3*82% 

670*62 

3.84% 

”ST 

4419.73 

36.03% 

4901.75 

•36.14% 

5383.76 

36.22% 

5865.78 

36.30% 

6347.80 

36.36% 

^PLANT 

54.24% 

55.22% 

56.01% 

56. 70; 

1 

z 

57. 

27% 



Qr “ 

REACTOR HEAT RATE 



w ** 

MHD 

MHD NET ELECTRIC POWER ; 

u 

MHD 


W 

GT 

GAS TURBINE ELECTRIC POWER : 

'^GT 


W » 

ST 

STEAM TURBINE ELECTRIC POWER: 

’^ST 

Si 


_ j ”mhd 

«GT 


^PLAHT 1 ^ 

Qr ■ 



^mhd " "^""compressor 

""gx ""compressor 

""sT ^ ""pump 


w. 


ST 


w. 


X 100 


MD 


W. 


100 + 


GT 

5r 


100 + 


ST 


100 
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Table 10.3 Plant Not Overall Bfflciencles For MUD Inlet Temperature Variation 


MHD Inlet 
Temperature 


3000°K 


3250”k 


3500“k 


4973.45 100.0% 5138.94 



319.12 6.42% I 319.12 6.2: 


1112.20 22.36% B 1112.20 21.64% 



3750°K 


5299.94 

100.00% 

2139.78 

40.37% 


5458.27 100.0% 



6.02% B 319.12 5.85% 


20.99% 


5693,55 

100*0% 

2590.04 

45.49% 

.319*12 

5.60% 

1112.20 

19.53% 


'PLANT 


62.75% 


65.11%! 


67.38% 


69.56% 


70.62% 


Qj^ “ REACTOR HEAT RATE 

Wmhd “ MHD NET ELECTRIC POWER; = «MHD " ^COMPRESSOR 

Wgj » GAS TURBINE ELECTRIC POWER; OUTPUT - 2 x 

Wgj •» STEAM TURBINE ELECTRIC POWER: Wg^ = Wg^ OUTPUT - Wp^jjjp- 


'plant 


+ ^ + X 100 

% ^R 


w \ / w u 

100 + “ 100 + 

^R \ \^R 




















Plant Net Ef f iciency (%) 



MHD Inlet Temperature 

1 Nuclear MHD Power Plant - 1 (with Heat Exchanger) 

2 Nuclear MHD Power Plant - 2 (with Mixing Tank) 

Fig. 10.6 Plant Net Efficiencies vs MHD Inlet Temperature 
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For Nuclear MHD Power Plant with regeneration, the major contribution 
of the electric power is produced in the top of the power cycle by the 
MHD subsystem (33*97% 45*49% from 100% heat produced by the reactor). 

The power production has been shifted toward the top of the ternary 
cycle with a positive effect on overall efficiency* This system produces 
overall efficiencies that are 25-35% higher than actual power plants in 
use and that are 15-20% higher than the expected coal^fired MHD power 
plants* 
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